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Many extensions of the Standard Model lead to non-standard neutrino interactions (NSI), which 
can affect the interpretation of current and future neutrino data. We explore an alternative ap- 
proach to the study of these dimension eight four fermion interactions (of coefficient 4eGf/V^), by 
bringing into play collider data. In models where coefficients of potential dimension six operators are 
suppressed by cancellations, LEP2 provides interesting bounds on NSI operators (e <, 10~ 2 — 10~ 3 ). 
If NSI are contact interactions at LHC energies, they induce an effective interaction qqW + W e+e^ 
(according to the Equivalence Theorem), to which the LHC at 14 TeV and with 100 fb _1 of data 
has excellent sensitivity (e^3x 10 -3 ). 
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I. INTRODUCTION 

Many extensions of the Standard Model, such as Su- 
persymmetry or leptoquarks, naturally induce low energy 
contact interactions of the form: 

e f a f^CfYPxf){v alpVf} ) , (1) 

where / € {u, d, e} is a first generation charged Standard 
Model fermion, and a, ft € {e,/i,r}. These are referred 
to as Non Standard neutrino Interactions (NSI) [1], and 
can be generated by SU(3) x SU{2) x U(l) gauge invari- 
ant operators at dimension six or higher. Future neutrino 
facilities, such as a Neutrino Factory [2], could be sensi- 
tive to such interactions with e <; 10 _4 [3-5]. Formalism 
and current bounds [6] on NSI are reviewed in section II. 

This paper aims to make a preliminary exploration of 
the complementarity of current collider experiments and 
future neutrino facilities to these neutral current opera- 
tors. We focus on neutral current NSI induced at dimen- 
sion eight (because the dimension six operators are more 
strictly constrained), from operators such as 

^{qYPLq){Hl alp H^) (2) 

where q and i are SM doublets, and H is the Higgs. 
Requiring e ^ 10~ 4 implies A 8 ^ 2 TeV. This suggests 
that the new mediating particles (Z's, leptoquarks, spar- 
ticles...), if weakly coupled or contributing in loops, are 
kinematically accessible to the LHC. In this case, their 
discovery prospects are model-dependent, and have been 
widely studied[7]. Here, to retain some degree of model 
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independence, we consider effective operators at collid- 
ers (LEP2, the Tevatron, the LHC). So for the LHC, 
we are assuming heavy New Physics with <; 1 couplings. 
We also assume that the New Physics does not gener- 
ate "dangerous" dimension six operators involving two 
charged leptons instead of two neutrinos, because these 
are strictly constrained. A catalogue of such models can 
be found in [8, 9]. 

In section III, we argue that in some models, such 
as those considered by Gavela et.al. [9], the dimen- 
sion eight NSI interaction of the form (1) is accompa- 
nied by dimension eight contact interactions involving 
charged leptons rather than neutrinos, with coefficients 
~ s/A$,(t — u)/A|, where s, t and u are the Man- 
delstam variables. We estimate that the LEP2 bounds 
on such contact interactions could be translated to e ^ 

icr 2 ,icr 3 . 

In section IV, we use the Equivalence Theorem to re- 
place (H)v a — > W + e a , and study the prospects for de- 
tecting qq — > W + W~£+£p at the LHC. Rough estimates 
suggest that couplings <; 1 are required for NSI to func- 
tion as dimension eight contact interactions at LHC en- 
ergies (assuming e <; 1CT 4 ). Such contact interactions 
would induce few events (<r(pp —> W = W~l + l~) <~ 1(V 3 
fb x(10~ 4 /e) 2 ), but at very high px where Standard 
Model backgrounds are negligeable. We estimate that 
the LHC could be sensitive to e <; 3 x 1(V 3 . 



II. REVIEW AND NOTATION 

Non Standard neutrino Interactions (NSI) are four- 
fcrmion interactions involving at least one neutrino, 
which are induced by Beyond the Standard Model 
physics. For a recent review, see e.g. [8]. In this paper, 
we focus on the neutral current NSI, sometimes called 
"matter NSI". At energies -C mw, the contact inter- 
actions between a neutrino, anti-neutrino and two first 



2 



Z couplings 


t f 

9l g R 


e, [i, t 
u, c, t 
d, s, b 


* , 

— \ + sin 2 9w sin 2 9w 
§ - § sin 2 9 W - 1 sin 2 9 W 
— \ + | sin 2 614/ | sin 2 9w 



TABLE I: Z couplings to SM fermions. 



generation fermions can be written 

Csm = -2V2G F 9 f P [^l P PL^}[h P Pf] (3) 
P,f,P 

Cnsi = -e f j2V2G F (9 alp P L iy p )(h p P x f) (4) 

where / is a first generation fermion, Px = Pl or Pr, 
and the tree-level g P are given in table I. We neglect SM 
loop corrections [10] in this paper; for a discussion which 
includes them, see [11]. 

NSI are phenomenologically interesting because they 
could interfere in the extraction of neutrino oscillation 
parameters [3, 4, 12, 13], clouding, for instance the mea- 
surement of sin6>i3 at a v Factory [14]. There have been 
several studies of the sensitivity of current [15] and fu- 
ture [5, 16] neutrino facilities to (charged current) NSI 
in the production and/or detection of neutrinos, and to 
neutral current NSI which could contribute a "matter" 
effect in neutrino propagation. Future facilities, such as a 
;/Factory, could be sensitive to e ^ 10 . Recent bounds 
on the coefficients of neutral and charged current NSI can 
be found in [6], and are of order 0(1) — 10~ 2 . 

From a more theoretical perspective, NSI are also in- 
teresting because they would indicate New Physics in the 
lepton sector at the TeV scale. However, dimension > 4 
operators induced by New Physics, such as Eqn. (4), 
should be invariant under the Standard Model gauge 
symmetries. It is not simple to construct New Physics 
models which would be discovered via NSI. For instance, 
the combination {p a l p vp) can easily be obtained as a 
component of {£ a j p £p)- However, an effective operator 
containing this current also gives rise to contact interac- 
tions involving the charged leptons (e a , y p ep), which are 
more strictly constrained(see e.g. [17]). Neutral current 
NSI, without undesirable charged lepton contact interac- 
tions, can be obtained at dimension six from operators 
such as 



{£H*)p{Hl 



(5) 



that modify the neutrino kinetic terms after spontaneous 
symmetry breaking [18]. However, such modifications to 
the normalization of the neutrino kinetic terms are con- 
strained by the unitarity of the PMNS matrix, implying 
that e 10~ 3 [18]. Alternatively, a dimension 6 contact 
interaction between neutrinos and electrons 



can be obtained by the tree-level exchange of a scalar 
with flavour antisymmetric couplings {{e.g. a slepton in 
R-parity violating Supersymmetry). Recent constraints 
on the coefficient of the operator (6) can be found on [8] . 

As noted in [19], gauge invariant operators that induce 
neutral current NSI, without accompanying charged lep- 
ton contact interactions, can be obtained at dimension 
eight, for instance via the operator 



TA {qYPLq){Hl alr Hl l3 ) 



(7) 



where q and £ are SM doublets, 
Requiring e Z 10 -4 implies 



At 



and H is the Higgs. 



(8) 



where v = (H) = 174 GeV, so A 8 <, 2 TeV. The obvious 
question arises: what New Physics could generate such 
a large coefficient at dimension 8 while remaining con- 
sistent with all other constraints? Various options have 
been explored. 

New Physics that induces the dimension 8 operators 
at tree level, without generating tree level dimension six 
operators, was explored in [8]. The models studied in 
[8] also induced charged current NSI, and non-unitarity, 
which restricted the es below their phenomenological 
bounds. However e ~ 10~ 3 could be obtained. A less 
restrictive approach to tree level New Physics for neutral 
current NSI was followed in [9] , where cancellations were 
allowed in the coefficient of the unwanted dimension six 
operators. For instance, in a model containing a singlet 
scalar and a singlet vector leptoquark, of masses ms and 
my and couplings to first generation fermions 



h(qHT 2 £)Sl + g(q^£)Vf 



(9) 



it is a straightforward exercise in Fierz transformations 
to show that the sum of the diagrams in figure 1 gives 



h 2 



2m\ 4m| 



[{qiq){h£) + {qifq){£in)) , (10) 



so parameters can be chosen such that the coefficient is 
zero. This cancellation will no longer be exact at non- 
zero momentum transfer. In section III of this paper, we 
show that LEP2 bounds on contact interactions generi- 
cally impose e ^ 10 -2 — 10 -3 on models which circumvent 
the dimension 6 bounds via such a cancellation. 

New Physics could also induce NSI via loop diagrams. 
This option is disfavoured by naive power counting: 



16tt 2 A4 



£ 

i)2 



(ii) 



V;](V*]) f 



(6) 



which implies that such New Physics should have a 
mass ~ 300 GeV - - accessible to current colliders, 
and constrained by precision electroweak data. How- 
ever, the attraction of generating neutral current NSI 
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at one loop, is that the SM has a built-in mechanism 
to suppress Flavour Changing Neutral Currents (FCNC) 
up to dimension eight: the (quadratic) GIM mecha- 
nism (Recall that GIM allows FCNC with coefficients 
oc m'jGp/(16Tr 2 ), where m/ = yjv is the mass of an in- 
ternal fermion, so one can interpret that the diagram has 
two additional Higgs legs.). Minimal Flavour Violation 
[20] extends the GIM mechanism to New Physics mod- 
els, thereby allowing NP with flavoured couplings at the 
TeV scale, despite the stringent bounds on quark FCNC. 
So one could anticipate that NSI will "naturally" arise 
via loops in models with TeV New Physics in the lepton 
sector. An obvious example is Supersymmetry, which 
can have spartner masses in the 300 GeV range while re- 
specting precision elextroweak constraints. However, su- 
persymmetric loop contributions to NSI were computed 
in [11], and it appears difficult to obtain e 10~ 4 . 



and where the dimension six four-charged-lepton oper- 
ator vanishes due to a cancellation (see e.g. figure 1, 
with the quarks replaced by electrons). Such models were 
classified in [9]. In sections IIIB and IIIC, we estimate 
bounds on A 8 . 
A. Dimension 8 operators with 4 charged leptons 

Gavela et al [9] classified the simplest extensions of the 
Standard Model which satisfied the following conditions: 

1. the New Physics (NP) generates neutral NSI at tree 
level 

2. the NP couples to SM bilinears (this allows to iden- 
tify the quantum numbers of the simplest combina- 
tion of NP that generates a dimension 6 or 8 oper- 
ator) 



U',D',Q 



U,D,Q 




U°,D°, Q 



U°,D°, Q 



FIG. 1: Scalar and vector leptoquarks can separately induce 
dimension six operators involving two quarks and two charged 
leptons, with a relative negative sign. The mass and coupling 
of the vector can be tuned relative to the scalar, to obtain a 
cancellation at zero momentum transfer, which allows these 
leptoquarks to participate in generating significant dimension 
8 NSI operators. 



III. LEP2 BOUNDS ON DIMENSION 8 
DERIVATIVE OPERATORS 

In this section, we estimate bounds from LEP2 on con- 
tact interactions involving four charged leptons, with co- 
efficients oc {s,t, w}/Ag. Some models which give rise to 
NSI also induce such interactions, so these bounds may 
be translated, in a model- dependent way, to the e coef- 
ficient of NSI involving two neutrinos and two electrons 
(/ = e in eqn 4). In section III A, we show that the 
dimension eight four-charged-lepton operators can arise 
in models where the NSI are generated at tree level, 



3. the four charged lepton operators must be sup- 
pressed, so either the NP does not induce them, 
or the various NP coefficients sum to zero 

As an example where the coefficient of a dangerous 
four-charged lepton operator vanishes due to a cancella- 
tion, consider a model containing an SU(2) doublet vec- 
tor of mass my and an SU(2) doublet scalar of mass 
with couplings 



h(et)Sl+g(e^t)Vf 



(12) 



Neglecting the masses of the external fermions, the dia- 
grams of figure 1 give the operators 



m v — u 



(e^rmh^ c ) 



2(m 2 s - t) 

g 2 
m v — u 



(efe)^) (13) 



where t and u are the Mandelstam variables. At zero 
momentum transfer, the coefficients of the dimension 6 
operator can be arranged to cancel by choosing g 2 /m v = 
h 2 /{2m 2 s ). However, at dimension 8, will arise operators 
such as 



g 2 {t/m 2 s -u/m v ) 



nit 



(14) 



Since s + 1 + u = (in the limit of neglecting the fermion 
masses), we take our double-derivative dimension eight 
operators to be oc s and cut — u: 



^((D'l^iD^^L^L) -+ JL(e a rP L e p )(e^P Y e) (15) 

8 8 

^((D^ a )7%)(^P^))-^ I (47' t (^)^(^(^i)) -+ ^(e^PLe^e^Pye) (16) 
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where L is a first generation lepton, either a doublet 
(Py = Pl) or singlet (Py = Pr), and to the right of 
the arrow is the four charged lepton component of the 
operator on the left. 



B. Bounds on flavour diagonal lepton operators 

LEP2 searched for dimension six contact interactions 
of the form 



±^{eYPxe)(f al ,P Y f a ) 

A 6,± 



(17) 



in the channels e+e~ — > e+e~, t+t~, qq (where 

q £ {it, d, s, c, 6} is a light quark; we postpone the quark 
case to section HID), by comparing the observed cross- 
section 

-l 



dcg 



da 
dee 



and forward-backward asymmetry 



Jo" 1 dn 



(Icq 



- 1+1 dc < 



da 



(18) 



(19) 



to SM expectations, eg = cos 9 is the angle between the 
incoming and outgoing fermions. 

To translate the published bounds on A 6j ± into con- 
straints on dimension eight operators cx s/A|, we focus 
on LEP2 data at large y/s : 183 -> 208 GeV. As there was 
significant initial state radiation, several experiments re- 
quired that the invariant mass of the final state leptons 
exceed 85% of the centre of mass beam energy. So we 
approximate 

(20) 



47T 



A 2 At Af 

where v = (H) = 174 GeV. If we assume that the NSI 
operator with coefficient i> 2 /A| is accompanied by the 
four charged lepton operator with coefficient s/A|, then 
we can extrapolate bounds of order 



e ^ 



471V 

iv 6 



3.8 x 10 



2 H 



(21) 



Consider now the dimension eight four charged lepton 
operator with coefficient (t — u)/A|, eqn (16). Since 



t — U = CgS 



(22) 



it is clear from eqn (19) that this operator will have dif- 
ferent contributions to a and Afb than those of eqns 
(15) and (17), because even powers of c 2 , in the matrix- 
element-squared contribute to cr, and odd powers to Afb- 
We focus on the interference term A4gM-M* A + h.c. 
between the SM and the contact interaction 1 , and we 



1 This term is generically larger than |.Ma| 2 , for A6 of order the 
experimental bounds. 





bound 


£ 


e + e~ —¥ e + e~ 
XY=LL 
LL 
RL 
RL 


A 6+ k 10.3 TeV 
A 6 _ k 8.3 TeV 
A 6+ k 8.8 TeV 

A 6 _ k 12.7 TeV 


<, 3.7 x 10~ 3 
<, 5.6 x 10~ 3 
<, 4.7 x 10~ 3 
<, 2.4 x 10 -3 


e + e _ — > 
XY=LL 
LL 
RL 


A 6+ k 8.1 TeV 
A 6 _ k 9.5 TeV 
A 6 ± k 6.3 TeV 


<, 5.9 x 10~ 3 
<, 4.3 x 10~ 3 
£ 9.1 x 10 -3 


e + e~ — > t + t~ 
XY=LL 
LL 
RL 
RL 


A 6+ k 7.9 TeV 
A 6 _ k 5.8 TeV 
A 6+ k 6.4 TeV 
A 6 _ k 4.6 TeV 


<, 6.2 x 10~ 3 
<, 1.1 x 10~ 2 
£ 9.1 x 10~ 3 
<, 1.8 x 10 -2 



TABLE II: The middle column gives the LEP2 bounds [21- 
23] on the dimension six contact interaction of the first col- 
umn, normalized according to eqn (17). These bounds can be 
used to constrain dimension eight operators with two deriva- 
tives and four charged leptons. If we assume that NP in- 
ducing the derivative operators also induces NSI, we obtain a 
bound on the e coefficient of the dimension eight NSI opera- 
tor, which is obtained from eqn (20), and given in the third 
column. Although the bounds are quoted with two signif- 
icant figures, they are merely order of magnitude estimates, 
as various constants could appear in the passage between four- 
charged-lepton-derivative operators and NSI operators. 



treat it as a small correction to the SM contribution. If 
we define the total cross-section in the presence of an 
operator of the type (15) to be a = <Jsm{^ + Acr s ), and 
Afb — Afb,sm(1 + AAp b,s), then some straightforward 
algebra 2 shows that 



Acr t _„ 



± Ao- s 
4 

±4AA 



FB,s 



(23) 
(24) 



where + and — are respectively for for Py — Pl , Pr in 
eqn (16). So to obtain bounds on four-charged-lepton- 
operators with coefficients cx (t — w)/A|, we should refit 
to the data. Instead, since the LEP2 data agrees with 
the SM, we assume that the bounds on the (t — «)/Af 



2 Recall that, for massless external fermions, t = — s(l — eg)/ 2 and 
u = — s(l + C(v)/2. To obtain a term in which contributes 

to Afb ( is linear in eg), from an operator of the form eqn (15) 
or (17), requires the contraction of a e a( g 7 5 from both the Dirac 
Traces on the initial and final state fermions. Since e a p e s arise 
from Tr{7„ 7^ 7575}, it is clear that a cx (gy + g A 2 )(gy + g A 2 ) 
and Afb °c 9y9 A 9y9 A , wnere t is the final state lepton. For 
operators cx (t — u) of the form of eqn (16), these relations are 
reversed: A FB cx (gy 2 + g A 2 )(g^ + 9 A ) and a « 9 V 9 A 9 V 9 A - 
Expanding the various g v A couplings as a SM part plus a NP 
part gives equations (24). 
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(n»P x e)(£^P Y l) 


bound 


e 


e + e" — > e ± fi T 
V XY 


A 8 £ 570 GeV 


& 8.7 x 10~ 3 


e + e~ — > e ± r T 
V XY 


A 8 £ 485 GeV 


<, 1.6 x 10~ 2 


e + e~ -¥ 
V XY 


A s £ 500 GeV 


<, 1.5 x 10~ 2 



TABLE III: Bounds on the scale of dimension eight four- 
charged-lepton derivative operators, normalized as eqns (15) 
and (16), from OPAL constraints [24] on lepton flavour- 
changing contact interactions (e + e — > e+e^ , a / P). These 

bounds come from |A1a| 2 , so are the same for all fermion 
chiralities and are independent of the sign of the contact in- 
teraction coefficient. The bounds on e given in the last colomn 
are estimated by assuming that the dimension eight NSI op- 
erator is induced with an 0(1) coefficient at the same scale. 



operators of cqn (16) are comparable to the bounds on 
the s/A| operators of cqn (15). 

Using bounds on k§ obtained by ALEPH [21] and 
OPAL [22] , and by Bourilkov from a combination of data 
from the four experiments [23], we use eqn (21) to esti- 
mate the bounds in table II. 



Flavour off-diagonal lepton operators 



The OPAL experiment saw one 



— > e ± ^ event 
at ^ = 189 - 209 GeV, and published limits [24] on 
a(e + e~ — > e ± n T ,e ± T T ,T ± [i T ) at LEP2 energies. These 
bounds give more stringent limits on operators of the 
form (15) and (16) than the LEP1 bounds on lepton 
flavour-changing Z decays. For y/s = 200 — 209 GeV, 
OPAL found 



cr(e+e -> e±ej) < < 



22 fb [e/x] 

78 fb [er] (25) 
{ 64 fb M 



We estimate that the cross-section induced by the oper- 
ators of eqns (15) and (16) is 




(26) 



and obtain the bounds on A 8 listed in table III (we 
neglect the 8/3, so estimate the bound from er 8 ~ 
s 3 / (327rA|).). If we assume that the four-charged-lepton- 
operators with coefficients oc s/Af , (t — u)/A|, would be 
accompanied by an NSI operator with coefficient u 2 /A|, 
then we can translate to bounds on e = u 4 /Ag, which are 
given in the last column of table III. 



(e^P L e)(q^P Y q) 


ALEPH bound on A 6 


e& 


Y=L 
Y=R 


A 6 £ 8.0+,9.7_ TeV 
A 6 ;> 5.2+, 4.1_ TeV 


<, 4.8 x 10~ 3 
& 1.8 x 10~ 2 



TABLE IV: The second colomn gives bounds from ALEPH 
[21] on the scale of dimension six contact interactions oper- 
ators, normalised as eqn (17). The estimated bounds on e 
given in the last colomn are obtained using eqns (20) and 
(21). 

D. Two-quark, two lepton operators 

LEP2 set bounds on dimension six contact interaction 
of the form of cqn(17), involving two electrons and two 
light quarks q G {u, d, s,c, b}. The bounds [21] were ob- 
tained assuming the same contact interaction to all quark 
flavours but the top. Nonetheless, we estimate that we 
can apply them separately to contact interactions involv- 
ing uorrf quarks. This is because the bounds arise mostly 
from the interference of the dimension six contact interac- 
tion with the SM matrix element, and the cancellations in 
this flavour sum make it of the same order for one or five 
flavours. The ALEPH [21] bounds on dimension six con- 
tact interactions are listed in the second colomn of table 
IV. We can translate these into bounds on the dimension 
eight double derivative operators of eqns (15,16) by ap- 
proximating s <~ v 2 ~ t — u, as discussed in section IIIB 
for lepton final states. We gucstimate the bounds using 
the average of A 6 +, A 6 ~, because the effect of the sign 
would be different for u and d type quarks. If we further 
assume that the dimension eight double-derivative oper- 
ators are accompagnied by the NSI operator with a co- 
efficient of the same order, then we can estimate bounds 
on e using eqn (20), which are given in the last colomn 
of table IV. 

Bounds on contact interaction involving two muons 
and two first generation quarks, normalised according to 
eqn (17), were obtained at the Tevatron by CDF [25], 
who searched for excesses in the lepton mass spectrum 
in 110 pb- 1 of data. Using 1 fb' 1 of data, DO [26] per- 
formed a "model-independent" search for new physics in 
a wide range of final states, including e ± T =F ,e ± /i =F and 
^ ± r =F , and did not find evidence for contact interac- 
tions. It is difficult to extrapolate either the CDF or 
DO bounds to dimension eight double-derivative contact 
interactions. CDF obtained bounds A 6 ^ 3 — 4 TeV 
for (uj a Pxu)(p,'y a PYn) interactions, and A 6 <; 1.5 — 2.3 
TeV for (drf° ' Pxd)(fir) a PY f-i) interactions. We calculated 
3 the contribution to <r([pp — > jj^i) arising in the presence 
of a contact interaction with coefficient s/Ag, and com- 
pared it to the contribution due to a contact interaction 
with a coefficient 47r/A§. Requiring that the dimension 
eight contribution be less than the allowed dimension six 



with CTEQ10 partem distribution functions 
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contribution gives Ag ^ 500(200) GeV, for contact inter- 
actions involving u(d) quarks. If one makes the further 
assumption that the dimension eight double-derivative- 
operators accompany NSI operators, this suggests that 
bounds of order e & 1 — .01 could be obtained for NSI 
involving first generation quarks and any leptons other 
than fr. 



IV. LHC DISCOVERY REACH 

In this section, we study the prospects at the LHC 
for NSI involving quarks, in the improbable, but com- 
partively mo del- independent, scenario that they appear 
as contact interactions. We first argue in section IV A 
that NSI can give a final state with W and charged 
leptons, then we estimate the production cross-section 
at the 14 TeV LHC, and finally discuss sensitivity and 
backgrounds in section IV C. 

One might "expect" the new particles which mediate 
NSI to be kinematically accessible to the LHC, so they 
would not induce contact interactions. Consider first NSI 
which arise at one loop, then e <~ w 4 /(167r 2 A 4 ) <; 10~ 4 
implies that the new particles of mass < A arc light 
enough to be produced at the LHC. Alternatively, if NSI 
arise at tree level, then e ~ v 4 /A 4 10~ 4 implies A ^ 2 
TeV. The new particles mediating NSI on quarks must 
couple to first generation quarks, so they must have a 
large mass, beyond the LHC reach, to avoid being di- 
rectly produced. The LHC reach for mediators such as 
Z's or leptoquarks, depends on the details of the partic- 
ular model, but is of order 3-5 TeV [7]. However, more 
realistically, the parameter 1/ A 4 is a product of coupling 
constants divided by a product of masses, so masses be- 
yond the LHC reach could arise for new couplings > 1, or 
if only some of the new particles are "sufficiently" heavy 
( e.g. A 4 = M 2 m 2 , M > m). In section IVD, we discuss 
a model with couplings < 1 and a hierarchy in masses, 
and show that the cross-section can be enhanced, be- 
cause some of the new particles can be produced on-shell 
so there is less phase space suppression. Despite the first 
sentence of this paragraph, we assume NSI parameters 
such that the only signal at the LHC would be a contact 
interaction. 



A. The Equivalence Theorem and NSI 



(E, 0, 0, k)/mw and longitudinal oc (k,0,0,E)/mw po- 
larisation states of the W ± gauge bosons. Since the time- 
like polarisation is "unphysical" , and cancelled by the 
goldstone, this suggests that the longitudinal can be 
replaced by the goldstone. 

More formally, after fixing the gauge of the electroweak 
sector, there are remaining global symmetries, which give 
BRS conditions, such as (in Feynman t'Hooft gauge) : 

d^W+v - m w H+ = . (27) 

These relations generate the electroweak Slavnov- Taylor 
Identities. We interpret the Equivalence Theorem to say 
[28] that this condition can be imposed on physical ma- 
trix elements. In particular, in an NSI operator, we would 
like to replace v a {H) — > Wlz&- 

As an example of how the Equivalence Theorem works, 
it is instructive to consider a one generation type I seesaw 
model, containing a heavy singlet majorana neutrino Nr, 
of mass Mr >• v, with a Yukawa interaction XNrH£, 
where A = 1. If SU(2) is unbroken, it is straightfor- 
ward to verify [29] that Nr has equal branching ratios 
to H°i>, H°*V, H + e~ and H~e + . After electroweak sym- 
metry breaking, the heavy mass eigenstate becomes 

L^Nr+^1 . (28) 

Despite being mostly singlet, with a very small admixture 
of i/£, L decays about half the time to a final state con- 
taining a W ± : BR(L -> W+e~) ~ BR(N R -> H+e~) = 
1/4 [30]. This arises because the l/m w in the W spin 
sum 

S"" = -g» v + ( 29 ) 
m w 

cancels the g 2 v 2 that appears upstairs in the squared 
matrix element: 

j(^) ^M.^}^ • (30) 

This confirms that at high energy, in the SSB theory, 
rates involving the W are similar to those for the gold- 
stone in the unbroken theory. We therefore assume that 
NSI operators induce qq — > W + W~e+e~p at the same 
rate as qq — > H + H~ e+e~p . 



In a gauge invariant dimension eight NSI operator, the 
HQH a v a Vfi interaction is accompagnied by H + H~e a ep, 
which could be expected to reincarnate, after electroweak 
symmetry breaking, as a vertex involving W + W~e a ep. 
This expectation can be formalised, at energies ^> m w , 
via the Equivalence Theorem[27, 28], which identifies the 
Goldstone with the longitudinal component of the 
W ± . 

Intuitively, the Equivalence Theorem relies on the 
similarity, at energies 3> rnw, of the timelike oc 



B. The production cross-section 

In the Equivalence Theorem limit, where the W ± is 
approximated as a scalar H ± , the matrix element for 
qq — > e+e~p H + H~ is comparatively simple, and the par- 
tonic cross-section is 



a(qq — > H + H~ 



n 



i i 

15 25 



A! 



2§ 2 



3(4tt^ 



(31) 
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This result was obtained by analytically integrating over 
massless four body phase space up to the last two inte- 
grals, which are performed by Maple™ [31]. It corre- 
sponds to 1/15 of a dimensional analysis estimate: the 
last parenthese is massless four-body phase space [32], 
and the square brackets enclose an estimate of the matrix 
element squared, together with a (27r) 4 which accompa- 
gnied the four-momentum (S-function. 

Integrating (31) over CTEQ10 parton distribution 
functions [33], assuming the contact interaction cou- 
ples equally to u and d quarks, we obtain <r(pp — > 
H+H-e+e~) at the 14 TeV LHC plotted in figure 2. The 
cross-section is small — suppressed by four body phase 
space — so the LHC could have difficulty detecting NSI 
that appear as contact interactions. 



g(W+W-l+l-) trom NSI at LHC|b ^ 14. r = 7), CTEQ10 (Q=100 GeV) | 




0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

A(TeV) 



FIG. 2: Estimated cross-section o(pp — > W + W e^e g ) in 
fb, at the LHC with 14 TeV (blue) and 7 TeV (red) centre of 
mass energy, due to dimension eight contact interactions with 
coefficient 1/A 4 . 

We obtain a(pp — > W + W~ e+e^) "by hand", via 
the Equivalence Theorem, because it is straightforward 
and manageable. We are unclear on the coefficient and 
structure of the correspondingg qqW + W~e+eZ interac- 
tion which would arise in the spontaneously broken elec- 
troweak theory, and which could be implemented in pro- 
grammes such as Madgraph [34]. The Madgraph cross- 
sections we obtain, using the interactions discussed in the 
following sections, are consistent 4 (to within a factor of 
a few) with figure 2. 



4 This agreement may be partly accidental. The "longitudi- 
nal W" component of the |matrix element| 2 for the vertex of 
Eqn. (32) is (p w + ■Pyy~) 2 /myy x the |matrix element| 2 for 
{cj'y^ q){e^^l tl/ ep)H Jr H~ . If Eqn. (27) implies that these |matrix 
element| 2 s should be equal, then the coefficient of Eqn. (32) has 
the curious value m 2 1/ /(A|p w ,+ • p w -). Using that the average 
(s) in the pp — > e^e^H^ H~ process is {§) ~ 18 TeV 2 , and that 

s = {P W + +Pw- +Pe+ +P e -) 2 . we § uess Pw+ ■ Pw- - L5 
TeV 2 . With this value, the Madgraph cross-section for the in- 
teraction of Eqn. (32) agrees with figure 2. 

In the model of section IV D, we also obtain a suppression of 
the qqW + W~ e& eg interaction with respect to the coefficient 



C. Detecting NSI contact interactions at the LHC 

To determine the reach of the LHC in e Q( g, we 
need to discuss backgrounds. The SM cross-section for 
W+W-£+£- is about 2 fb at ^fs = 14 TeV. Besides 
irreducible backgrounds, one should also consider back- 
grounds faking our signal. For example, as discussed in 
[35], W ± Z/"/* and jets has a cross-section of 30 fb at 
the 14 TeV LHC, and could look like W + W~£+£p when 
the jets can be reconstructed to a W ± . More overwelm- 
ing could be tt+jets, which could produce signals such as 
trileptons, where both W's decay leptonically and one of 
the &-jets fakes an isolated lepton. To study the effect of 
a 6-jet faking a lepton we took a 3M tt sample generated 
with ALPGEN [36] using MLM matching [37] and then 
showered with PYTHIAv6.4 [38] at 14 TeV. The cross 
sections for tt + (0,1,2) jets are (440,778,730) pb. The 
efficiency of asking for 3 leptons in the sample is about 
5xl0 -3 , which is consistent with the estimates in [39]. 
This effect would reduce the tt background to levels of 
100's of fb, still too large for most of the NSI parameter 
space (see figure 2). 

Fortunately, asking for multi-leptons is not our only 
handle on reducing backgrounds. We anticipate that an 
effective operator will produce highly boosted objects, 
allowing an efficient cut on px, and also that the NSI 
events may be spherical. To study these expectations, 
we simulate the simple interaction 

{qi P P L q){W+T alp W->n+) (32) 

using FeynCalc [40] and MadGraph v5 [34]. q can be 
a u or d quark, and both are assumed to couple with 
the same strength. We impose parton level cuts of 1 77 1 < 
2.5 and px > 20 GeV. If we attribute to this interaction 
the coefficient of the interaction (35), with Mr = Ag, 
MadGraph gives the same total cross-section as section 
IV B. 

In Fig. 3 we plot the px of objects in both the ti and 
signal samples. The signal and background distributions 
are well separated, and asking for a cut on px of both lep- 
tons of order 400 GeV would reduce the ti backgrounds 
to less than 10 -5 fb and keeps 70% of our signal. We can 
interpret that the s 3 dependence of eqn (31) encourages 
energetic final states. Indeed, integrating s under the 
cross-section gives an average ^/s for NSI events of order 
four TeV, or one TeV per particle, in agreement with the 
Pt distribution obtained in the simulation. 



of the gauge invariant NSI operator. This Im^lM?^ factor has 
a fixed mass in the denominator (this is reassuring). However, 
Mr is not determined by Ag, and we do not expect to need ad- 
ditional details of the model to obtain cr(qq — > e^e^W+IV") in 
the Equivalence Theorem limit. Setting Mjj = A ~ TeV, gives a 
coefficient numerically similar to the curious estimate above for 
the coefficient of Eqn (32), but this could be an accident. Notice 
that we did not compute the | matrix element | 2 for Eqn. (32). 



The pt cut is very effective for removing the SM back- 
ground. So one could relax the assumption on the num- 
ber of electrons and muons, and search for NSI contact 
interactions by requiring high px multiparticle events. In 
particular, this would give sensitivity to NSI operators in- 
volving taus. However, if such events were seen, boosted 
taus would be hard to tag, although they may show up 
as fat jets. 




P T (GeV) 

FIG. 3: Distribution of pr of the leptons in NSI (blue) and ti 
(black). The normalization of the distributions is chosen to 
show both signal and background in the same scale, and the 
energy of the LHC is 14 TeV. Note that the vertical axis is 
logarithmic. 



is characterized by well-separated, highly boosted objects 
in a pretty spherical event. Using these characteristics, 
especially a cut on p T , we showed that the largest back- 
ground, from ti can be reduced below the signal. The 
LHC reach depends on the value of e, but if we assume 
the NSI signal is background-free, and ask for 100 events 
at a luminosity C (in fb^ 1 ), then the reach in e is 

£ - 3 x 10- 2 /y/C (33) 

For example, for a luminosity of 100 fb _1 , the 14 TeV 
LHC could be sensitive to NSI-induccd contact interac- 
tions corresponding to 

e > 3 x 10~ 3 . (34) 

This estimate can be compared with the e > 10 sensi- 
tivity of a neutrino factory. 

So far we have discussed the LHC run at 14 TeV, but 
one can already obtain bounds using the current run at 
7 TeV. The reduction in cross section from 14 to 7 TeV 
leads to a decrease in cross section of about a factor 60. 
Therefore, assuming the current running would end up 
with 30 ftr 1 of combined ATLAS+CMS data, the sensi- 
tivity with the 2011-12 data set would be e ~ 4 x 10~ 2 . 

D. A Z' model 



In Fig. 4 we show the sphericity of the signal versus the 
ti sample. As expected, the signal is more spherical than 
the background. The signal objects are well separated in 
AR space, so no issues of overlapping between the leptons 
and the W's arise. The signal is not particularly forward: 
the effect of \rn\ < 4.5 to 2.5 is reducing the signal by less 
than 8%. 
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FIG. 4: Sphericity in the signal (blue) and ti (black). The 
normalization of the distributions is chosen to show both sig- 
nal and background in the same scale. 



In this section, we discuss how an interaction similar to 
Eq. 32 could arise in a toy model for NSI, and compare 
this Z' model predictions with the study done in the 
previous section. 

Let us assume a new Z', with coupling g' q = 1 to first 
generation quarks, so it is produced at hadronic colliders 
via qq — > Z' . We suppose a cancellation at the level 
of dimension-6 operators (as in Gavela et al. [9]), so 
neglect the Z' contribution to the dijet cross section. Our 
Z' also couples, with g' N = 1, to new singlet leptons 
Nr , of mass Mr, which have a yukawa coupling A to 
SM leptons 5 . After electroweak symmetry breaking, the 
heavy mass eigenstate L will contain a small admixture of 
doublet neutrino (see eqn (28)), and will have a coupling 
cx Au/Mr to Ws. This allows the process illustrated in 
Fig. 5. In the contact interaction limit for the Z' and L 
propagators, this diagram gives the interaction 



9 ? 

9 v 
A 4 M 



2 {qi P P L q){t a WlpWh) 



(35) 



which gives a somewhat more complicated matrix ele- 
ment than the interaction simulated in the previous sec- 
tion (see eqn (32)). 

We simulated the process of Fig. 5 with MadGraph, 
requiring a large mass for the Z' resonace, to ensure that 



Summarizing, the cross-section is small, but the signal 



It is unlikely that this toy model is consistent with low energy 
neutrino data. 
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FIG. 5: The high energy origin of pp — > e%e^ W + W in a Z' 
model with new leptons. 

it is produced off-shell. The L mass was allowed to vary 
in the v — ¥ TeV range. This allows to explore an example 
where some of the new particles responsable for NSI (the 
Nr), are kinematically accessible to the LHC. The final 
state kinematics for this process are very similar to those 
of the effective operator analysis in the previous section. 
In Fig. 6 we show the px distribution and sphericity for 
the effective operator and a model with a 3 TeV Z' and 
lighter N4 with mass of 300 GeV. 

Although the kinematic distributions are very similar 
in the Z' model and the contact interaction model, the 
cross sections are dramatically different. Indeed, after 
integration over PDFs, the cross sections are roughly re- 
lated by the two-body versus four-body phase factors [32], 

o z , ~ 10 5 x a s (36) 

which would lead to a larger NSI signal. If we ask for 100 
events at a luminosity L (in fb _1 ), then the reach in e is 

£ ~ 10~ 4 / \TL (37) 

which for a luminosity of 100 fb _1 would lead to a sensi- 
tivity in the range of e <; 10~ 5 . 




50 if 
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FIG. 6: Comparison of the kinematic distributions for the Z' 
model and the effective operator. 



V. SUMMARY/CONCLUSION 

In this paper, we studied some implications at colliders 
of neutral current, Non-Standard neutrino Interactions 
(NSI) of the form of eqn (1). We assumed these con- 
tact interactions to be the remnants, after elextroweak 
symmetry breaking, of gauge-invariant dimension eight 
operators such as eqn (2). Notation and models that 
could give such interactions are reviewed in section II. 

In section III, we considered the class of models where 

• these operators are induced at tree level by New 
Physics 

• this New Physics does not induce "dangerous" di- 
mension six operators involving charged leptons, 
because the various contributions to their coeffi- 
cients cancel. 

We argued that this cancellation only obtains at zero mo- 
mentum transfer, so the New Physics could induce con- 
tact interactions involving charged leptons rather than 



neutrinos, with coefficients ~ s/A| (where s in the cen- 
tre of mass energy). These contact interactions must 
satisfy LEP2 (and Tevatron/LHC) bounds, which can 
be translated via eqn (20) into approximate bounds e 
10" 2 , 10~ 3 on the coefficient of NSI.JThe bounds on NSI 
operators involving (ee)(£ a £ a ), (ee)(£ a £b) for a ^ /3, and 
(ee)(qq), are given respectively in tables II, III, and IV. 
It is possible to constrain flavour-changing four-lepton in- 
teractions because OPAL published bounds [24] on con- 
tact interactions giving e ± fi^ z , t fi^, and e ± r =F in the 
final state. 

In section IV, we focussed on NSI that would behave 
like a contact interaction at LHC energies. We used the 
Equivalence Theorem to relate NSI operators involving 
Higgses and neutrinos, with operators with two charged 
gauge bosons and two charged leptons. The cross-section 
we obtain for pp — > W + W~ e+e^ is small (see figure 2), 
due to four-body final state phase space suppression, but 
gives spherical events with highly-boosted objects. On 
the other hand, section IV D discusses a Z' model where 
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some of the new particles can be produced on-shell, which 
enhances the the cross section by orders of magnitude. 
Interestingly, the final state kinematics remain similar. 

We showed in section IV C that asking for multilep- 
tons, and in particular, imposing a hard cut on pr (see 
figure 3), can reduce Standard Model backgrounds below 
the signal. This is even the case for NSI involving rs. 
Our results suggest that the LHC with 100 fb _1 of data, 
could be sensitive to contact interactions induced by NSI 
with e^3x 10~ 3 (e <; 10~ 5 ) for the model independent 
(dependent) case, (see eqns (33) and (37)). 
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